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ABSTRACT 


One -dimensional  laser -Doppler  velocimetry  measurements 
were  taken  with  standard  optics  in  back  scatter  mode  across  a 
normal  shock  at  a  Mach  number  of  1.35.  Back  pressure  on  a 
blow-down  supersonic  tunnel  was  controlled  to  place  a  normal 
shock  in  a  4  by  4  inch  test  section  and  schlieren 
visualization  techniques  were  used  to  verify  and  record  shock 
position  and  behavior.  Velocity  surveys  were  taken  across  the 
shock,  using  various  filtering  techniques,  in  an  attempt  to 
quantify  shock  unsteadiness.  Additional  surveys  were  performed 
to  further  characterize  the  flow  in  the  test  section.  The 
velocity  surveys  upstream  and  downstream  of  the  shock  compared 
favorably  with  pressure  and  temperature  data  and  normal  shock 
relations.  Surveys  across  the  shock  indicated  distinct  and 
repeatable  velocity  patterns,  and  the  measured  location  of  the 
shock  matched  schlieren  photographs. 
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I.  INTRODUCTION 


A.  THEORETICAL  BACKGROUND 

In  many  experimental  studies  in  fluid  mechanics  and 
aerodynamics,  accurate  velocity  measurements  within  flowfields 
are  of  fundamental  importance.  Detailed  research  of  complex, 
turbulent,  or  high  speed  flows  is  virtually  impossible  without 
the  aid  of  sophisticated  velocity  measuring  devices.  Laser 
Doppler  velocimetry  (LDV)  is  a  non- intrusive  optical  technique 
that  uses  the  Doppler  principle  to  measure  the  velocity  of 
particles  in  a  fluid.  Over  the  past  three  decades,  LDV  has 
evolved  into  a  highly  adaptive  and  precise  tool  that  is  widely 
used  in  the  measurement  of  flow  fields. 

The  relative  motion  between  a  source  of  radiation  and  a 
moving  body  will  cause  the  observed  frequency  of  the  source  to 
be  Doppler  shifted.  When  a  beam  is  directed  onto  a  particle 
moving  in  the  flow,  the  particle  will  scatter  light.  The 
scattered  light  is  frequency  shifted  and  the  velocity  can  be 
measured  from  the  shift  in  this  scattered  light. 

Although  there  are  several  types  of  optical  systems 
utilized  by  LDV,  the  most  common  type  is  the  dual  beam  system 
which  uses  optics  to  divide  the  laser  beam  into  two  parallel 
beams  of  equal  intensity.  The  two  beams  are  then  focused  by  a 
lens  such  that  they  intersect  at  a  point  in  the  flow.  The 
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constructive  interference  from  the  coherent  light  waves  of 
each  beam  form  interference  fringes  and  as  a  particle  passes 
through  the  fringes  it  will  alternately  scatter  or  not  scatter 
light  as  it  crosses  a  light  or  dark  fringe. 

The  rate  at  which  the  particle  crosses  the  fringes  is 
equal  to  the  rate  of  oscillation  of  the  signal  from  the 
scattered  light  which  is  collected  by  a  lens  and  directed  onto 
a  photo  detector.  Fringe  spacing  can  be  calculated  by  the 
equation: 


S- 


X _ 

2  sinG 


(1) 


where  X  is  the  wavelength  of  light  used  and  6  is  the  half 
angle  of  the  two  beams  after  passing  through  the  focusing  lens 
[Ref.  1]  .  By  obtaining  the  frequency  of  the  oscillating 
signal,  fd,  with  a  photo  detector,  the  velocity  of  the 
particle  can  be  determined  by  the  following  equation  [Ref.  1]  . 


\fD 

V  =  Sfn  =  - r— 77 

°  2  sin0 


(2) 


If  the  particle  is  small  enough  such  that  it  moves  with  the 
flow  of  the  fluid  then  the  velocity  of  the  fluid  is  equal  to 
the  velocity  of  the  particle. 

Frequency  shifting  is  commonly  employed  as  a  means  of 
creating  a  fringe  velocity.  A  Bragg  cell  is  used  to  shift  one 
of  the  dual  beams  thus  creating  a  pattern  of  moving  fringes. 
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The  fringe  movement,  with  the  appropriate  downshifting  of  the 
signal,  allows  measurement  of  reverse,  low  velocity,  and 
highly  turbulent  flows. 

The  region  of  intersection  of  the  two  beams  is  referred  to 
as  the  probe  volume.  Each  beam  must  have  a  Gaussian  intensity 
distribution  and  the  volume  created  when  the  two  beams 
intersect  is  an  ellipsoid.  As  a  particle  crosses  the  probe 
volume  two  separate  signals  are  produced.  The  Doppler  signal, 
as  discussed  earlier,  is  produced  by  the  relative  motion  of 
the  particle  with  respect  to  the  fringes.  Also,  a  pedestal 
signal  is  created  due  to  the  Gaussian  variation  of  the  light 
intensity  in  the  beams. 

The  Doppler  signal  is  superimposed  on  the  pedestal  signal 
which  can  easily  be  filtered  out  since  the  pedestal  frequency 
ranges  from  0  to  several  kHz  and  the  Doppler  frequency  is  in 
the  range  of  hundreds  of  kHz  to  Mhz .  The  total  signal  that  is 
produced  by  a  particle  crossing  the  probe  volume  is  termed  a 
Doppler  burst.  The  duration  of  the  burst  can  be  estimated  by 
knowing  the  velocity  of  the  particle  and  the  dimensions  of  the 
probe  volume.  It  is  these  bursts  that  are  processed  to  obtain 
velocity  information.  [Ref.  1] 

The  signal  to  noise  ratio  (SNR)  of  a  LDV  signal  is  given 

by: 

SNR  =  [  D*D*'2  ]  2dD2GV2  (3) 

A  f  iaf  p 
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where  r?q  is  the  quantum  efficiency  of  the  photo  detector,  PQ 
is  the  laser  power,  At  is  the  velocity  bandwidth,  dp  is  the 
particle  size,  G  is  the  scattering  parameter,  V  is  the 
visibility,  Da  is  the  collection  aperture,  ra  is  the  focal 
length,  De_2  is  the  diameter  of  the  intersecting  region,  and 
C  is  a  constant  of  proportionality  [Ref.  2] .  It  can  be  seen 
from  the  above  equation  that  laser  power  or  purticle  size  can 
be  increased  to  increase  the  SNR  but  power  increases  can  be 
expensive  and  larger  particles  may  not  follow  the  flow 
accurately.  Increasing  t ne  collection  aperture  by  increasing 
the  size  of  the  collecting  lens  or  increasing  the  diameter  of 
the  intersecting  region  by  means  of  beam  expansion  are  more 
popular  methods  of  increasing  the  SNR. 

The  beam  that  is  emitted  from  the  laser  source  is 
vertically  polarized  and  has  a  divergence  angle  of  a  few 
milliradians .  It  is  essentially  a  parallel  beam  but  also  has 
a  region  of  minimum  diameter  which  is  referred  to  as  the 
waist.  The  light  waves  at  the  waist  are  plane  waves  therefore 
it  is  here  that  beam  crossing  is  accomplished.  If  chis  is  not 
the  case,  the  fringes  that  are  formed  are  curved  instead  of 
parallel  which  can  cause  errors  in  velocity  measurement.  To 
correct  for  this  phenomenon,  most  systems  install  collimators 
before  the  beams  are  split  to  increase  the  parallelism  of  the 
beams . 

The  source  of  the  LDV  signal  is  totally  dependent  on  the 
particles  that  are  entrained  in  the  flow.  Although  particles 
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that  are  naturally  present  in  most  flows  are  capable  of 
producing  LDV  signals,  artificial  particle  generation  is 
generally  used  to  produce  particles  of  uniform  size.  It  is 
highly  desirable  to  have  particles  which  have  good  reflective 
properties  and  the  ability  to  closely  follow  the  flow.  Size 
becomes  a  critical  factor  in  that  the  particles  must  be  large 
enough  to  scatter  light  effectively  yet  small  enough  to  avoid 
lagging  through  areas  of  high  velocity  gradients. 

The  introduction  of  particles  to  a  flow  is  termed  seeding 
and  is  accomplished  by  atomizers  and  fluidized  bed  generators. 
There  are  numerous  types  of  seeding  material  which  can  be  used 
depending  on  the  application.  Atomized  water  or  oils  can  be 
used  when  uniform  particle  size  is  not  critical.  Solid 
particles  in  solution,  such  as  polystyrene  latex  (PSL) 
particles  in  an  alcohol  solution,  can  be  utilized  to  create 
particles  of  a  more  uniform  size.  This  type  of  solution  is 
introduced  at  a  location  such  that  the  alcohol  evaporates  in 
the  flow  prior  to  the  test  section  leaving  only  the  solid 
particles . 

Photo  detectors  convert  the  scattered  light  into  a 
voltage.  This  voltage,  which  corresponds  to  the  Doppler 
frequency,  is  used  by  the  signal  processor  to  obtain  the 
velocity  data.  One  of  the  more  common  types  of  signal 
processors  used  in  LDV  systems  is  the  counter.  A  counter 
measures  the  frequency  of  the  Doppler  burst  by  measuring  the 
time  required  for  a  number  of  cycles.  For  example,  the  time  is 
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measured  for  8  cycles  in  a  signal.  This  time  is  then  compared 
with  the  time  for  5  cycles  of  the  same  signal.  The  ratio  of 
the  two  times,  in  an  ideal  case,  would  be  5  to  8.  An  error 
limit  is  set  on  the  ratio  as  a  means  of  validation  and  if  the 
time  is  within  the  preset  limits,  it  is  used  to  calculate  the 
velocity.  [Ref.  2] 

The  obvious  advantage  that  LDV  presents  over  conventional 
techniques  such  as  pitot  probes  and  hot  wires  is  the  ability 
to  make  measurements  without  disturbing  the  flow  field.  LDV  is 
also  ideally  suited  for  precise  measurement  of  single  or 
multiple  components  of  the  velocity  vector,  needs  no 
calibration,  and  can  be  used  in  a  wide  variety  of  fluids  and 
flows.  Despite  many  advantages,  however,  LDV  does  possess  a 
unique  set  of  problems  such  as  the  dependence  on  effective 
seeding,  the  cost  of  optical  systems  and  signal  processors, 
and  the  difficulty  of  correctly  interpreting  LDV  data  by 
statistical  analysis. 


B .  HISTORICAL  BACKGROUND 

LDV  was  first  used  to  measure  the  velocity  of  a  flow  field 
by  Yeh  and  Cummins  in  1964  [Ref.  3].  Since  that  time,  a 
continuous  effort  has  been  made  by  the  Aerodynamic  community 
to  develop  and  improve  LDV  as  an  invaluable  means  of  accurate 
flow  velocity  measurement.  As  Stevenson  [Ref.  4]  points  out 
in  his  historical  review  of  LDV,  advances  have  been  focused  on 
optical  design,  signal  processing,  and  data  analysis. 
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The  first  system  used  by  Yeh  and  Cummins  split  the  light 
emitted  from  a  helium- neon  laser  into  a  measurement  and  a 
reference  beam,  as  shown  in  Figure  1.  This  configuration  was 
similar  to  a  Mach-Zehnder  interferometer  but  was  difficult  to 
use  due  to  the  inability  to  consistently  align  the  various 
optical  components.  Around  1970,  systems  began  to  incorporate 
a  differential  Doppler  arrangement  which  divided  the  beam  into 
two  parallel  beams  that  could  be  focused  at  a  point  by  a 
single  lens  as  shown  in  Figure  2.  Hanson  [Ref.  5]  theorized 
that  a  variation  in  the  fringe  spacing,  which  would  effect  the 
Doppler  frequency,  was  a  potential  problem  due  to  the  Gaussian 
nature  of  the  laser  beam.  Durst  and  Stevenson  [Ref.  6] 
verified  Hanson's  theory  experimentally  and  eliminated  the 
unwanted  variations  by  placing  correction  lenses  between  the 
laser  and  the  beam  divider. 

Another  optical  advance  has  been  the  practice  of 
decreasing  the  focused  beam  diameter  by  means  of  expanders 
inserted  between  the  beam  divider  and  the  focusing  lens  as 
shown  in  Figure  3.  Also,  many  advances  have  been  made  in  the 
general  construction  of  precise  optical  equipment  which  has 
led  to  refined  LDV  measurements. 

Improvements  in  signal  processing  units  have  also  played 
a  large  role  in  the  development  of  LDV  systems.  Devices  such 
as  the  frequency  tracker,  burst  processor  (counter) ,  and  the 
photo  detector  have  been  specifically  designed  to  handle  the 
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signals  encountered  in  LDV  measurements.  Early  work  with 
frequency  trackers  proved  reliable  up  to  about  50  MHz  and  with 
heavy  seeding.  Modern  counters,  however,  allow  the  user  to 
select  the  number  of  cycles  to  be  counted  and  have  a  Doppler 
frequency  limit  of  approximately  200  MHz.  Photo  detectors 
were  first  introduced  as  a  means  of  obtaining  measurements 
from  low  levels  of  scattered  light  around  1972.  Improvements 
have  continuously  been  made  to  the  frequency  capabilities  and 
processing  time  of  these  units.  [Ref.  7] 

Data  analysis  can  be  divided  into  two  general  categories. 
First,  corrections  for  velocity  errors  associated  with  seeding 
and  probe  volume  size  effects.  Second,  corrections  for  large 
errors  associated  with  turbulence  due  to  processor 
characteristics  and  the  nature  of  the  Doppler  signal  [Ref.  4]  . 
Methods  for  correcting  the  velocity  errors  have  included 
averaging  to  obtain  a  mean  particle  velocity  and  rejecting 
high  amplitude  signals  to  eliminate  large  particle  lag  time. 
McLaughlin  and  Tiederman  [Ref.  8]  used  time  averaging  rather 
than  particle  averaging  to  eliminate  velocity  error.  They 
theorized  that  more  high  velocity  than  low  velocity  particles 
would  cross  a  probe  volume  over  a  given  period  of  time  which 
would  lead  to  an  incorrect  mean  velocity  reading.  Methods  used 
to  correct  for  errors  in  measured  turbulence  parameters  have 
become  well  understood  and  are  corrected  by  methods  described 
by  George  [Ref.  9]  . 
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Refinements  over  the  last  twenty  years  have  allowed  LDV  to 
become  an  integral  part  of  a  wide  range  of  aerodynamic 
studies.  LDV  measurements,  such  as  those  performed  by  Goebel, 
Dutton,  Krier,  and  Renie  [Ref.  10],  have  documented  the 
behavior  of  supersonic  mixing  layers .  Supersonic  separated 
flow  in  a  compression  corner  has  been  investigated  with  the 
use  of  three-dimensional  traversing  LDV  systems  as  in  the  work 
of  Baroth  [Ref  11] .  Detailed  studies  of  flows  through 
turbomachinery  have  been  performed  such  as  the  tests  conducted 
by  Ceman  [Ref.  12]  in  which  he  used  LDV  to  study  air  flows 
through  a  high-speed  ducted  fan.  Precise  airfoil  evaluations 
and  studies  of  trailing  edge  flows  have  also  been  presented  as 
in  the  report  of  Absil  and  Passchier  [Ref.  13].  Virtually 
every  facet  of  aerodynamic  flow,  at  velocities  ranging  from 
near  stationary  to  well  into  the  supersonic  range,  has  been 
examined  with  the  aid  of  LDV. 

One  area  that  has  been  studied  intensely  over  the  past 
five  years  is  shock- boundary  layer  interaction.  The  current 
designs  for  turbofan  engines,  which  are  called  upon  to  deliver 
high  levels  of  thrust,  require  fan  and  leading  compressor 
stage  relative  Mach  numbers  that  are  supersonic.  The  location 
and  control  of  the  shock  systems  in  such  stages  is  of  great 
interest  to  the  designers.  A  shock  will  likely  form  at  the 
leading  edge  of  each  blade  and  will  impinge  on  the  suction 
side  boundary  layer  of  the  adjacent  blade  as  shown  in  Figure 
4.  This  shock  structure  is  termed  a  lambda- foot. 
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Figure  4.  Shock- Boundary  Layer  Interaction,  From  Ref.  18 
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[Ref.  14]  In  order  to  understand  and  better  control  the 
losses  associated  with  this  type  of  shock  system,  detailed 
velocity  measurements  must  be  obtained  in  a  non- intrusive 
manner.  LDV  is  perfectly  suited  for  this  type  of  measurement. 

Another  aspect  of  LDV  that  has  been  studied  in  detail  is 
particle  dynamics  in  the  flow  field.  Since  the  basis  of  LDV 
techniques  is  the  measurement  of  particle  velocity,  it  is 
essential  that  accurate  models  are  developed  to  predict 
particle  behavior  in  regions  of  instability  or  high  velocity 
gradients.  Studies  done  on  particle  lag  prediction,  Chesnakas 
and  Andrew  [Ref.  15],  and  particle  dynamics  effects  on  LDV 
measurements,  Bloomberg,  Dutton,  and  Addy  [Ref.  16],  have 
attempted  to  verify  the  accuracy  of  theoretical  models  used  to 
predict  particle  behavior  through  planar  oblique  shocks. 

C .  PURPOSE 

The  purpose  of  the  present  work  was  to  establish  the 
capability  to  make  LDV  measurements  in  transonic  flow.  This 
would  allow  additional  studies  to  be  conducted  in  the  area  of 
shock -boundary  layer  interaction.  The  objectives  of  the 
present  work  focused  on  three  areas.  First,  the  existing 
supersonic  tunnel  located  in  the  Gas  Dynamics  Laboratory  (GDL) 
of  the  Naval  Postgraduate  School  was  refurbished, 
instrumented,  and  modified  to  allow  introduction  of  seeding 
material.  Secondly,  a  traversing  LDV  system  capable  of 
accurate  one -dimensional  velocity  measurements  was  set  up. 
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Finally,  experimental  LDV  measurements  were  taken  in  the  free 
stream,  through  the  test  section  boundary  layer,  and  across  a 
normal  shock. 
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II.  EXPERIMENTAL  APPARATUS 


A.  SUPERSONIC  WIND  TUNNEL 

The  wind  tunnel  used  was  a  blow- down  type  supersonic 
tunnel  located  in  the  GDL  at  the  Naval  Postgraduate  School .  A 
photograph  and  schematic  of  the  tunnel  are  shown  in  Figures  5 
and  6,  respectively.  Supply  air  was  produced  by  a  3  stage 
centrifugal  compressor,  dried,  and  stored  at  150  psi  in  an 
8000  cubic  foot  tank  farm.  A  supply  air  schematic  is  shown  in 
Figure  7. 

Supply  air  at  150  psi  was  controlled  to  the  wind  tunnel  by 
means  of  a  pneumatic  control  valve.  Additional  supply  air  was 
regulated  at  12-14  psi  and  routed  to  the  valve  to  allow  manual 
control .  A  pressure  probe  and  absolute  pressure  gauge  were 
installed  at  the  top  center  of  the  plenum  to  allow  precise 
control  of  plenum  pressure.  The  plenum  consisted  of  a 
cylindrical  chamber  approximately  20  inches  in  diameter 
containing  a  9  inch  diameter  circular  flat  plate  mounted 
perpendicular  to  the  flow  as  shown  in  Figure  8.  The  plenum 
contained  no  other  screens  or  flow  straightening  devices. 

The  contraction  and  circular  to  rectangular  transition 
section  had  an  area  ratio  of  5.7:1.  Sets  of  interchangeable 
aluminum  blocks  could  be  attached  downstream  of  the 
contraction  to  form  the  convergent -divergent  nozzle  and  test 
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CONVERGENT  -  DIVERGENT 
NOZZLE 


EXHAUST 


Figure  6.  Schematic  of  Supersonic  Wind  Tunnel 
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Figure  7.  Supply  Air  Schematic 


Figure  8.  Plenum  Conf igurat ion 
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section.  The  blocks,  as  shown  in  Figure  9,  were  comprised  of 
an  identical  pair  which  were  held  in  place,  top  and  bottom,  by- 
two  flat  plates  which  formed  the  side  walls  of  the  nozzle  and 
test  section.  There  were  four  sets  of  blocks  available  in  the 
laboratory,  each  designed  to  produce  a  specific  Mach  number  in 
the  test  section.  The  nominal  Mach  1.4  blocks  were  utilized 
throughout  this  study. 

The  test  section  measured  4  inches  horizontally  and  4 
inches  vertically.  Optical  access  was  provided  by  6  inch 
diameter  circular  glass  windows  on  each  side  of  the  test 
section  side  walls.  Pressure  ports  were  located  upstream  and 
downstream  of  the  test  section  and  metal  blanks  could  be 
installed  in  the  window  frames  to  provide  additional  pressure 
ports  across  the  test  section.  Pressure  transducers  were 
installed  just  upstream  of  the  test  section  and  on  the  plenum 
to  record  the  pressure  ratio  during  each  run.  Instrumentation 
and  software  for  these  transducers  is  discussed  in  Appendix  P. 
Air  was  exhausted  downstream  of  the  test  section  through 
ducting  to  atmosphere  outside  the  laboratory. 

Run  times,  based  on  constant  plenum  pressure,  varied  from 
a  minimum  of  2  minutes  with  plenum  pressure  maintained  at  30 
psig,  to  a  maximum  of  4  minutes  with  plenum  pressure 
maintained  at  15  psig.  Approximately  15  minutes  were  required 
between  each  run  to  allow  supply  air  to  build  up  to  150  psi. 

Ports  were  opened  downstream  of  the  test  section  to 
relieve  back  pressure  and  allow  the  shock  to  travel  past  the 


Figure  9.  Mach  1.4  Nozzle  Block 
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test  section.  A  valve  was  installed  on  one  of  the  ports  as 
ameans  of  controlling  the  back  pressure  and  positioning  the 
normal  shock. 

B.  LASER  DOPPLER  VELOCIMETRY  SYSTEM 

1.  Laser  and  Optics 

The  laser  and  optics  used  were  the  Lexel  model  95 
four-Watt  argon- ion  laser  and  the  TSI  four-beam,  two  color 
LDV  system.  The  system  was  mounted  on  a  TSI  aluminum  base 
which  was  attached  to  the  bed  of  a  commercial  milling  machine 
modified  to  serve  as  a  traverse  mechanism.  The  milling  machine 
could  be  manually  traversed  in  the  X  (streamwise) ,  Y 
(vertical) ,  and  Z  (horizontal)  directions  to  allow  surveying 
within  the  test  section.  A  photograph  of  the  laser,  optics, 
and  milling  machine  is  shown  in  Figure  10. 

Output  from  the  laser  was  initially  passed  through  a 
beam  collimator  to  ensure  that  the  waist,  (minimum  diameter  of 
the  beam)  ,  and  the  focal  point  coincided.  The  beam  then 
entered  the  color  separator  which  consisted  of  an  attenuator, 
dispersion  prism,  and  mirror  set.  The  attenuator  provided 
control  of  the  beam  intensity  and  the  dispersion  prism 
separated  the  laser  beam  into  the  two  strongest  olor  lines, 
the  514.5  nm  (green)  and  488  nm  (blue)  .  The  ~-i..  tor  set  then 
reflected  the  beams  out  of  the  color  separator  ;  ailel  to  the 
base.  For  the  present  study,  only  one -dimensional  velocity 
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Figure  10.  Laser,  Optics,  and  Milling  Machine 
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measurements  were  made  therefore  the  green  beam  was  blocked  at 
this  point  allowing  only  the  blue  beam  to  continue. 

After  leaving  the  color  separator  the  blue  beam  passed 
through  a  polarization  rotator  which  ensured  vertical  beam 
polarity.  The  beam  was  then  split  in  the  horizontal  plane  to 
allow  velocity  measurements  in  the  flow  direction.  A  single 
beam  was  frequency  shifted  with  the  use  of  a  Bragg  cell 
causing  the  fringe  pattern  to  move  in  the  direction  of  the 
flow.  The  frequency  was  shifted  40  Mhz,  but  downshifting  was 
not  used.  All  data  from  histograms  and  statistics,  therefore, 
had  to  be  corrected  for  the  40  Mhz  shift.  Both  beams  passed 
through  a  divergence  section  (beam  expander)  which  decreased 
the  probe  volume  and  increased  the  SNR.  An  end  lens  produced 
a  focal  length  of  762  mm.  [Ref.  17] 

The  optic  configuration  produced  a  3.1  degree  half 
angle  and  a  probe  volume  that  was  133  fim  in  diameter  and  2.5 
mm  in  length.  Fringe  spacing  was  4.51  /zm  with  28  fringes  in 
the  probe  volume.  A  TSI  model  9160  photomultiplier  system  was 
used  to  collect  the  backscattered  light  after  it  passed  back 
through  the  end  lens  and  beam  expander.  [Ref.  17]  A  schematic 
of  the  laser  and  optics  is  shown  in  Figure  11. 

2.  Data  Acquisition 

Data  acquisition  was  accomplished  with  a  TSI  model 
1990C  counter- type  signal  processor  unit.  The  unit  included  an 
input  conditioner,  timer,  digital  readout,  and  data  interface. 
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Signals  from  Che  photodeteccor  were  sent  to  the  processorwhich 
transformed  the  signal  into  voltages  proportional  to  the 
Doppler  frequency.  The  signal  processor  was  interfaced  with 
a  386  personal  computer  which  utilized  "FIND"  software 
provided  by  TSI  to  analyze,  present,  and  store  the  acquired 
data.  The  LDV  software  is  discussed  in  detail  in  Appendix  B. 
A  photograph  of  the  data  acquisition  system  is  shown  in  Figure 
12  . 


3 .  Seeding 

Seeding  was  accomplished  with  a  TSI  model  9306  six- jet 
atomizer.  Designed  specifically  for  LDV  applications,  the  TSI 
atomizer  consisted  of  a  liquid  reservoir,  pressure  regulator, 
atomizers,  dilution  system,  and  aerosol  outlet.  Regulated  air 
was  used  to  draw  liquid  droplets  from  the  reservoir  and 
impinge  them  on  a  spherical  impactor  which  caused  the 
formation  of  an  aerosol  that  exited  through  the  outlet.  [Ref. 
18] .  A  photograph  and  schematic  of  the  atomizer  are  shown  in 
Figures  13  and  14,  respectively. 

The  TSI  atomizer  was  capable  of  generating  particles 
from  virtually  any  liquid  or  suspended  particle  solution.  The 
majority  of  work  performed  in  the  present  study  was  done  with 
a  2%  solution  of  1  n m  sized  polystyrene  latex  (PSL)  particles 
suspended  in  alcohol.  The  particles  were  introduced  into  the 
flow  at  the  contraction  section,  as  shown  in  Figure  13.  The 
outlet  of  the  atomizer  was  tapered  and  fit  with  a  3/8  inch 
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Figure  12.  LDV  Data  Acquisition  System 


Figure  13.  TSI  Six- Jet  Atomizer 
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Figure  14.  Atomizer  Schematic,  From  Ref.  18 
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inside  diameter  aluminum  tube,  open  at  the  end,  and  extending 
into  the  tunnel.  The  open  end  of  the  tube  was  positioned  at 
the  center  of  the  flow.  Typically,  the  atomizer  was  operated 
at  a  supply  pressure  of  65  psi  with  0%  dilution  air. 

Care  had  to  be  taken  to  ensure  that  the  atomizer  was  not 
over  pressurized  either  from  the  supply  or  outlet  side.  Each 
situation  potentially  caused  spillage  of  seed  material  and  in 
the  latter  case,  prevented  seed  from  entering  the  contraction 
section.  To  avoid  both  situations  the  supply  pressure  was 
regulated  at  a  minimum  of  40  psi  (to  overcome  the  tunnel  back 
pressure  on  the  atomizer)  and  a  maximum  of  65  psi. 


C .  SCHLIEREN  SYSTEM 

A  schlieren  system  was  utilized  to  record  shock  position 
and  structure  within  the  test  section.  A  continuous  or  spark 
light  source  was  available  from  a  combination  unit.  The  light 
source  was  directed  to  a  parabolic  mirror  which  reflected  the 
beam  through  the  test  section.  The  light  was  then  passed 
through  a  collimating  lens,  filter,  and  knife  edge  arrangement 
before  entering  the  camera.  Most  photographs  were  impact 
shadowgraphs  taken  using  the  spark  source  without  utilizing 
the  filter  or  knife  edge.  A  diagram  of  the  arrangement  is 
shown  in  Figure  15. 
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III.  EXPERIMENTAL  PROCEDURE 


A.  OVERVIEW 

The  objective  of  the  present  study  was  to  perform  LDV 
measurements,  in  back  scatter  mode,  in  transonic  flow.  The 
experimental  procedure  was  divided  into  the  following  phases: 
apparatus  set-up,  schlieren  visualization,  and  LDV 
measurements.  Due  to  the  limited  run  time  of  the  supersonic 
tunnel,  set-up  procedures  were  critical  in  assuring  that  each 
run  was  valid,  well  documented,  and  efficient,  schlieren 
visualization  was  used  extensively  as  a  means  of  determining 
experimental  conditions  prior  to  LDV  measurements.  Three  types 
of  LDV  measurements  were  performed:  free  stream,  boundary 
layer  surveys,  and  normal  shock  surveys. 

One  of  the  main  objectives  during  the  experiment  was 
consistent  test  section  conditions.  After  the  initial  set-up, 
special  attention  was  given  to  operating  the  tunnel  under 
identical  conditions  for  each  subsequent  run.  This  was 
accomplished  by  keeping  plenum  pressure  and  back  pressure 
valve  position  constant  for  each  run  after  the  set-up. 
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B.  SET-UP  PROCEDURES 


1 .  Wind  Tunnel 

Two  set-up  runs  were  required  before  LDV  measurements 
could  be  performed.  The  tunnel  was  first  run  at  a  plenum 
pressure  of  15  psig  with  the  back  pressure  ports  full  open. 
This  produced  conditions  such  that  the  normal  shock  would 
travel  downstream  of  the  test  section.  The  passage  of  the 
normal  shock  was  observed  visually.  As  the  shock  passed  the 
test  section,  the  back  pressure  valve  was  manually  adjusted 
closed,  moving  the  shock  upstream,  until  the  shock  was 
positioned  in  the  test  section. 

Knowing  that  the  adjustments  to  the  back  pressure 
valve  would  restrict  the  flow  while  starting  the  tunnel  and 
cause  th^  stabilized  shock  position  to  change  (locate  further 
upstream)  when  the  tunnel  was  restarted,  the  shock  was 
visually  positioned  slightly  downstream  of  the  desired 
location.  The  tunnel  was  then  shut  down  and  run  again  with  no 
adjustments  to  the  back  pressure  valve  to  determine  the  shock 
position  for  subsequent  runs. 

Pressure  transducers  on  the  plenum  (P2)  and  test 
section  (Px)  were  calibrated  after  the  set-up  runs  using  the 
HP  Data  Acquisition/Control  system  and  "SPEED_5"  software  as 
described  in  Appendix  A.  Px  and  P2  were  then  calibrated  just 
prior  to  each  additional  run.  Several  runs  were  performed  with 
a  temperature  transducer  installed  in  the  plenum  to  measure 
the  total  temperature  which  was  used  to  calculate  static 
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temperature  in  the  test  section.  Test  section  static 
temperature  was  used  to  calculate  velocity  and  Doppler 
frequency  which  provided  information  used  to  adjust  filter 
settings  on  the  counter. 

2.  Schlieren 

The  set  up  of  the  schlieren  system  involved  checking 
the  alignment  and  spark  source.  Also,  a  wire  and  a  metal 
pointer  were  attached  to  the  exterior  of  the  test  section 
window.  Since  the  position  of  the  wire  and  pointer  and  the 
distance  that  separated  them  were  known,  the  position  of  the 
normal  shock  could  be  exactly  determined  in  all  photographs. 
A  small  free  jet  nozzle  was  used  to  focus  and  adjust  the 
schlieren  optics  prior  to  tunnel  operation. 

3.  LDV  and  Data  Acquisition 

Initial  alignments  were  performed  on  the  LDV  optics  to 
ensure  proper  beam  polarization  and  crossing.  Additional 
periodic  alignment  checks  were  performed  to  maintain  initial 
alignments.  Proper  beam  alignment  to  the  test  section  was  also 
checked  during  the  set-up  phase  to  ensure  that  the  beams  were 
oriented  parallel  to  the  flow  and  the  LDV  optics  were 
perpendicular  to  the  test  section. 

Data  acquisition  set-up  was  accomplished  utilizing 
"FIND"  software  as  described  in  Appendix  B  and  by  manually 
setting  processor  filters.  Cycles  per  burst,  and  timer 
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comparison  were  fixed  throughout  the  study.  The  LDV  and  data 
acquisition  systems  were  then  checked  by  operating  the 
atomizer  "open  air"  and  directing  the  probe  volume  into  the 
flow.  Data  rate  and  histogram  appearance  were  checked  to 
ensure  proper  operation  of  the  optics,  processor,  and 
software . 

As  previously  mentioned,  40  Mhz  shifting  was  used 
without  downshifting.  The  40  Mhz  shift  corresponded  to  a 
velocity  shift  of  180  m/s  with  fringe  spacing  4.512  /zm. 
Therefore,  all  velocity  data  obtained  from  the  LDV  data 
acquisition  system  was  corrected  by  adding  180  m/s.  Since  the 
values  for  turbulence  intensity  were  a  function  of  the  mean 
velocity,  turbulence  intensity  data  was  also  corrected. 

C.  SCHLIEREN  VISUALIZATION 

After  completion  of  set-up  procedures,  the  tunnel  was  run 
for  the  purpose  of  recording  shock  position  and  movement  by 
schlieren  visualization.  The  LDV  optics  were  removed  from  the 
field  of  view  of  the  schlieren  system  by  traversing  the  laser 
and  optics  base  below  the  test  section.  The  spark  light  source 
was  powered  and  tested  and  high-speed  film  was  preloaded  into 
the  camera. 

The  tunnel  was  allowed  to  start  and  stabilize  at  a  plenum 
pressure  of  15  psig  before  schlieren  photographs  were  taken. 
Multiple  photographs  were  shot  during  the  run  until  plenum 
pressure  dropped  below  15  psig.  The  visual  data  was  then 
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examined  and  recorded  as  a  means  of  verifying  shock  position 
and  determining  the  extent  of  shock  unsteadiness  (movement) . 

D .  LDV  MEASUREMENTS 

1.  Free  Stream 

Initial  LDV  measurements  were  made  upstream  and 
downstream  of  the  normal  shock  at  a  midstream  location.  The 
tunnel  was  started  and  "SPEED_5"  was  run  to  monitor  the 
plenum- to- test  section  pressure  ratio.  Seeding  was  introduced 
as  the  plenum  pressure  stabilized  at  15  psig.  All  six  atomizer 
jets  were  used  with  dilution  air  set  at  zero  and  supply 
pressure  set  at  65  psig. 

Several  single  point  measurements  were  taken  upstream 
to  verify  supersonic  flow  at  approximately  Mach  1.4  and 
several  single  point  measurements  were  made  downstream  to 
verify  subsonic  flow  at  approximately  Mach  0.8.  Processor 
filters  were  set  at  1 00  Mhz  (high)  and  20  Mhz  (low)  for 
upstream  locations  and  50  Mhz  (high)  and  10  Mhz  (low)  for 
downstream  locations.  Cycles-per-burst  was  fixed  at  4  and 
timer  comparison  was  fixed  at  5%  throughout  the  entire  study. 
Laser  power  was  also  a  constant  throughout  the  entire  study  at 
3  Watts. 

2 .  Boundary  Layer  Survey 

After  the  LDV  system  was  verified  as  accurate  during 
the  free- stream  measurements,  a  boundary  layer  survey  was 


performed.  The  tunnel  was,  again,  operated  at  15  psig  plenum 
pressure  but  a  change  was  made  to  the  seeding  location.  The 
outlet  tube  of  the  atomizer  was  lowered  in  the  contraction  to 
approximately  the  height  of  the  test  section  lower  wall  in  an 
attempt  to  introduce  more  seeding  particles  in  the  lower  half 
of  the  test  section. 

The  probe  volume  was  manually  positioned  as  close  as 
possible  to  the  lower  wall  of  the  test  section  at  a  point  well 
upstream  of  the  normal  shock  location.  After  the  tunnel  was 
started,  seeding  was  introduced  and  the  probe  volume  was 
traversed  manually  upwards  in  0.02  inch  increments  through  the 
boundary  layer.  Processor  filters  were  set  at  100  Mhz  (high) 
and  10  Mhz  (low)  during  the  entire  run.  Plenum  to  test  section 
pressure  ratios  were  monitored  throughout  the  run  and  the 
survey  was  ended  when  plenum  pressure  dropped  below  15  psig. 

3 .  Shock  Survey 

The  atomizer  outlet  was  positioned  in  the  middle  of 
the  contraction  section  and  a  shock  survey  was  performed  at  a 
midstream  location.  Operation  of  the  tunnel  was  identical  to 
the  boundary  layer  survey.  Seeding  was  introduced  after  the 
plenum  pressure  was  stable  at  15  psig  and  processor  filters 
were  set  at  100  Mhz  (high)  and  10  Mhz  (low)  .  The  survey  was 
started  at  a  point  upstream  of  the  normal  shock  and  traversed 
downstream  of  the  shock  in  increments  of  0.05  inches.  The  run 
was  stopped  when  plenum  pressure  dropped  below  15  psig. 
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IV.  RESULTS  AND  DISCUSSION 


A.  OVERVIEW 

The  results  of  the  present  study  are  divided  into  four 
areas:  schlieren  visualization,  free  stream  LDV  measurements, 
LDV  boundary  layer  surveys,  and  LDV  surveys  across  the  normal 
shock.  Selected  data  from  each  area  is  presented  as  it  applies 
to  the  overall  results  of  the  study.  Tabulated  data  appears  in 
Appendix  D. 

Wind  tunnel  runs  were  recorded  by  dating  each  run  and 
identifying  a  particular  run  on  a  particular  day  by  letter. 
For  example,  030293a  indicates  the  run  was  performed  on  March 
3rd,  1993  and  it  was  the  first  run  of  the  day.  All  data, 
(schlieren  photographs,  pressure  measurements,  LDV  files)  ,  are 
identified  with  a  particular  tunnel  run. 

Horizontal  position  in  the  test  section  was  defined  as 
inches  downstream  from  the  maximum  upstream  position  of  the 
probe  volume.  The  maximum  upstream  position  of  the  probe 
volume  was  restricted  to  1.5  inches  from  the  upstream  side  of 
the  test  section  due  to  traverse  table  limitations.  Vertical 
position  in  the  test  section  was  defined  as  inches  from  the 
bottom  of  the  test  section. 
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B .  SCHLIEREN  VISUALIZATION 


Initially,  schlieren  visualization,  using  a  continuous 
light  source,  was  utilized  to  view  the  normal  shock  movement 
during  tunnel  operation.  This  established  the  fact  that  the 
normal  shock  was  positioned  on  the  far  upstream  side  of  the 
test  section  with  plenum  pressure  set  at  30  psig.  Ports  were 
opened  downstream  of  the  test  section  to  relieve  the  back 
pressure.  With  the  ports  full  open,  the  shock  was  observed  to 
move  downstream  of  the  test  section. 

A  valve  was  installed  to  control  the  back  pressure  and 
video  recordings  of  the  schlieren  images  revealed  that  the 
position  of  the  shock  could  be  controlled  by  means  of  back 
pressure.  Closing  the  valve  (increasing  the  back  pressure) 
moved  the  normal  shock  upstream.  It  was  eventually  concluded 
from  multiple  test  runs  that  the  shock  could  be  positioned  in 
the  test  section  by  means  of  back  pressure  adjustment  with  the 
plenum  pressure  set  at  15  psig.  All  subsequent  runs  were 
performed  at  15  psig  plenum  pressure  which  approximately 
doubled  the  run  times. 

Spark  source  schlieren  photography  was  used  to  obtain 
information  on  the  movement  and  characteristics  of  the  normal 
shock  within  the  test  section.  Multiple  photographs  were  taken 
during  single  runs  to  establish  a  maximum  and  minimum  on  the 
extent  of  shock  movement.  Figures  16  and  17  show  10  spark 
source  photographs  taken  during  run  030393b. 
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The  photographs  clearly  identify  the  normal  shock.  Position 
and  scale  could  be  obtained  from  the  wire  and  pointer  which 
were  also  clearly  visible  in  the  prints.  The  wire  was  at  0.75 
inches  and  the  pointer  was  located  l  inch  downstream  from  the 
wire.  The  extent  of  the  normal  shock  movement  was  therefore 
determined  to  be  0.5  inches  with  a  maximum  upstream  position 
of  0.75  inches  and  a  maximum  downstream  position  of  1.25 
inches.  Although  the  position  of  the  shock  varied  from  run  to 
run,  the  extent  of  movement  was  found  to  be  almost  constant. 

Multiple  exposure  photography  was  attempted  during  several 
runs.  Figure  18  shows  a  triple  exposure  taken  during  030593a. 
An  upstream  position  of  the  shock  is  clearly  identifiable.  The 
downstream  shock  image  appears  darker  on  the  actual  photograph 
which  indicates  that  two  exposures  captured  the  shock  in 
essentially  the  same  position.  It  was  found,  however,  that  the 
print  became  overexposed  with  more  than  two  exposures , 
therefore  this  type  of  visualization  proved  to  be  of  limited 
value . 

The  resolution  in  most  of  the  photographs  was  adequate 
enough  to  accurately  determine  position  and  structure  of  the 
shock.  The  lambda  foot  of  the  normal  shock  was  clearly 
identified  in  some  photographs.  Figure  19  shows  a  well  defined 
lambda  foot  at  the  top  of  the  test  section  where  the  normal 
shock  interacted  with  the  boundary  layer. 
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Figure  18.  Triple  Exposure  Spark  Source  Schlieren 
Photograph,  Run  030593a 


Figure  19.  Spark  Source  Schlieren,  Run  030493a 


40 


C.  FREE  STREAM  LDV  MEASUREMENTS 


Several  upstream  and  downstream  LDV  measurements  were 
taken  at  a  midstream  location  with  the  shock  positioned  in  the 
middle  of  the  test  section.  These  runs  were  performed  to 
determine  the  accuracy  of  the  LDV  data  and  as  a  means  of 
finetuning  the  LDV  and  data  acquisition  systems.  Filter 
settings  on  the  processor  were  100  Mhz  (high)  and  20  Mhz  (low) 
upstream  of  the  shock  and  50  Mhz  (high)  and  10  Mhz  (low) 
downstream  of  the  shock. 

The  ratio  of  total  plenum  pressure  to  static  test  section 
pressure  (upstream  of  the  shock)  was  acquired  by  the  tunnel 
data  acquisition  system  and  processed  as  described  in  Appendix 
A  to  determine  upstream  Mach  number.  A  printout  of  the  tunnel 
data  acquisition  system  results  are  shown  in  Figure  20.  The 
data  indicates  a  rise  in  P2  (plenum  pressure)  as  the  tunnel 
was  started.  A  P2  reading  of  34  psia  corresponded  to  a  plenum 
pressure  of  20  psi  (as  read  from  the  plenum  pressure  gauge) . 
When  the  normal  shock  passed  the  upstream  side  of  the  test 
section,  as  indicated  by  the  abrupt  decrease  in  P1;  the  flow 
in  the  upstream  side  of  the  test  section  became  supersonic,  as 
shown  in  the  Mach  column  of  the  printout. 

Pressure  measurements  show  the  upstream  Mach  number  to  be 
stable  at  approximately  1.35.  Each  line  of  data  on  the 
printout  was  generated  in  about  5  seconds.  By  counting  the 
number  of  lines  in  the  printout  at  a  constant  plenum  pressure, 
a  rough  approximation  of  run  time  can  be  calculated. 
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Free- stream  LDV  measurements  indicated  an  upstream  velocity  of 
390  m/s.  By  knowing  the  value  of  plenum  stagnation  temperature 
(typically  48  F)  and  substituting  the  measured  LDV  velocity 
over  the  speed  of  sound  for  the  Mach  number,  the  test  section 
upstream  static  temperature  and  Mach  number  could  be 
calculated  as  follows: 


M  = 


V 


\jyRT 


(4) 


is  =  i  +  rs±M2 


(5) 


Tn  =  T  + 


V * 
2  Cr 


(6) 


T  =  Tn  - 


V * 
2  Cr 


(7) 


M  = 


V 


y 2 

yR(T0  -  JL-  ) 


2  C-r 


(8) 


The  Mach  number  calculated  from  LDV  measurements  and  the 
methodology  described  above  was  1.36  which  was  within  1%  of 
the  Mach  number  calculated  from  pressure  measurements. 
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Free -stream  LDV  measurements  indicated  a  downstream 


velocity  of  270  m/s.  Using  the  same  method  as  shown  above,  a 
Mach  number  of  0.86  was  calculated  downstream  of  the  normal 
shock.  The  expected  downstream  value  obtained  from  normal 
shock  tables  was  Mach  0.76.  The  difference  between  the 
measured  and  the  theoretical  values  is  at  least  partially  due 
to  the  boundary  layer  growth  in  the  test  section  downstream  of 
the  shock.  The  tunnel  was  poorly  designed  in  that  no 
allowances  were  made  for  boundary  layers  which  form  on  all 
walls.  The  boundary  layers  create  a  contraction  in  the  test 
section  which  must  cause  the  subsonic  flow  downstream  of  the 
shock  to  accelerate.  Because  of  this  effect,  a  Mach  number 
higher  than  that  shown  in  the  normal  shock  tables  could  be 
expected. 


D.  LDV  BOUNDARY  LAYER  SURVEYS 

Boundary  layer  surveys  were  performed  to  characterize  the 
flow  in  the  test  section.  All  surveys  were  done  between 
midstream  and  the  bottom  surface  of  the  test  section  with  the 
normal  shock  positioned  well  downstream.  The  probe  volume  was 
positioned  as  close  to  the  wall  as  possible  and  traversed 
upwards.  Overlapping  surveys  were  done  starting  below 
midstream  and  traversing  down  into  the  boundary  layer. 

Several  problems  were  encountered  while  conducting  the 
boundary  layer  surveys.  First,  limited  run  time  prevented 
taking  surveys  of  more  than  approximately  12  points .  This 
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limited  the  amount  of  survey  travel  with  small  increments 
between  each  point.  Secondly,  positioning  the  probe  volume 
:.ear  the  wall  by  visual  methods  was  less  than  exact.  Due  to 
the  nature  of  the  beams  and  the  theoretical  assumption  that 
the  velocity  measurement  is  taken  at  the  center  of  a  probe 
volume  with  finite  diameter,  obtaining  a  velocity  measurement 
at  the  wall  is  impossible. 

The  probe  volume  diameter  created  by  using  the  TSI  model 
9169-750  lens  was  133/nm  [Ref.  19].  If  the  probe  volume  was 
resting  on  the  wall,  the  velocity  measurement  were  actually 
taken  at  a  height  of  66.5/xm  or  0.0026  inches  (1/2  the  probe 
volume  diameter) .  Also,  since  the  beams  were  waisted,  the 
entire  probe  volume  had  to  be  positioned  at  a  finite  distance 
above  the  center  of  the  wall  in  order  for  the  beams  to  clear 
the  wall  as  they  entered  the  test  section.  This  distance  was 
approximately  0.0004  inches  [Ref.  19].  Thus,  when  the  probe 
volume  was  placed  as  low  as  possible  in  the  test  section,  the 
approximate  distance  from  the  wall  to  the  center  of  the  probe 
volume  was  .003  inches. 

Figures  21  and  22  show  the  results  of  boundary  layer 
surveys  031693a  and  031693b,  respectively.  In  031693a,  the 
traverse  was  upwards  from  the  wall  and  in  031693b  the  traverse 
was  downward  from  0.2  inches,  to  overlap  measurements  taken  in 
031693a.  Both  surveys  were  made  with  the  normal  shock 
positioned  downstream  of  the  test  section. 
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Figure  21.  Boundary  Layer  Survey  031693a 
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The  extent  of  the  boundary  layer  was  determined  by 
analyzing  not  only  velocity  distribution  but  also  turbulence 
intensity  changes.  The  boundary  layer  edge  velocity  was 
established  in  the  free  stream  LDV  measurements  as 
approximately  390  m/s.  An  examination  of  the  velocity 
distribution  in  031693a  and  031693b  indicated  a  boundary  layer 
thickness  of  between  0.18  inches.  The  turbulence  intensity 
measurements  taken  in  031693a  and  031693b,  which  should 
theoretically  increase  in  the  boundary  layer,  supported  the 
velocity  data  in  that  they  indicated  a  similar  boundary  layer 
thickness.  In  general,  the  data  obtained  from  the  boundary 
layer  surveys  was  repeatable  and  indicated  a  surprisingly 
small  turbulent  boundary  layer. 

The  data  obtained  in  both  surveys  were  normalized  and 
compared  with  theoretical  predictions  as  shown  in  Figures  23 
through  26.  Experimental  u/Ue  was  calculated  as  the  ratio  of 
measured  velocity  to  free-stream  velocity.  The  outer  edge  of 
the  boundary  layer,  <5,  was  defined  as  the  value  of  y  for  which 
u/Ue  was  0.99.  The  theoretical  value  of  u/Ue  was  calculated 
from  the  1/7  th  power  velocity  distribution  law  [Ref.  20]: 

l 

JL  =  r Z)  7  (9) 

Ue  5 

which  is  used  to  model  the  velocity  profile  in  the  boundary 
layer  on  a  flat  plate.  The  displacement  thickness,  6*,  which 
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Figure  23.  Normalized  Boundary  Layer  Data  and 

Calculations,  031693a,  Near  Wall  to  Free-Stream 
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Figure  24.  Boundary  Layer  Survey,  031693a,  Near  Wall  to 
Free-Stream 
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Figure  25.  Normalized  Boundary  Layer  Data  and 

Calculations,  031693b,  Free-stream  to  Wall 
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Figure  26.  Boundary  Layer  Sur 
Wall 
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represents  the  extent  of  the  displacement  of  the  external 
streamlines  due  to  the  growth  of  the  boundary  layer,  was 
calculated  by  the  following  equation  [Ref.  21]: 


The  momentum  thickness,  9,  which  is  an  alternate  definition  of 
boundary  layer  thickness  defined  in  relation  to  the  momentum 
deficit  in  the  boundary  layer,  was  calculated  from  the 
following  equation  [Ref.  21]: 

(11) 

The  form  factor,  H,  which  measures  the  fullness  of  the 
velocity  profile'  was  determined  by  taking  the  ratio  of 
displacement  thickness  to  momentum  thickness. 

In  both  surveys,  the  experimental  curves  agreed  with 
curves  generated  from  the  l/7th  law  model.  Experimental  data 
from  031693b  indicated  a  displacement  thickness  of  0.09  and  a 
momentum  thickness  of  0.07  which  yielded  a  shape  factor  of 
1.34.  Experimental  data  from  031693a,  however,  indicated  a 
displacement  thickness  of  0.11  and  a  momentum  thickness  of 
0.04  which  yielded  a  shape  factor  of  2.76.  Since  the  flow  was 
turbulent,  and  a  shape  factor  above  2.4  generally  indicates 
laminar  flow,  the  data  from  031693a  must  be  reanalyzed. 
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The  difference  between  the  two  surveys  was  that  in 
031693a,  a  velocity  measurement  close  to  the  wall  was  not 
obtained.  During  all  boundary  layer  surveys,  the  experimental 
velocities  were  biased  high.  By  adding  a  velocity  measurement 
at  the  wall,  biasing  it  low  (by  setting  it  equal  to  the 
theoretical  value) ,  and  recalculating,  the  results  were  in 
much  better  agreement  with  theory.  The  recalculations  and 
redefined  curve  are  shown  in  Figures  27  and  28.  Displacement 
thickness  was  recalculated  as  .08  and  momentum  thickness  was 
recalculated  as  .05  which  yielded  a  shape  factor  of  1.5. 

E.  LDV  NORMAL  SHOCK  SURVEYS 

Multiple  surveys  across  the  normal  shock  were  conducted  at 
a  midstream  location.  The  shock  was  positioned  in  the  test 
section  and  various  increments  were  used  to  traverse  upstream 
and  downstream  of  the  shock.  The  purposes  of  the  surveys  were 
to  determine  if  shock  position  and  movement  could  be 
accurately  predicted  with  LDV  measurements  and  to  study 
particle  dynamics  through  the  shock,  schlieren  visualization 
was  also  used  to  determine  shock  position. 

The  tabulated  results  from  surveys  030393a  and  031093b  are 
shown  in  Figures  29  and  30.  The  plotted  data  for  the  surveys 
appear  in  Figures  31  and  32.  In  both  surveys,  velocity  and 
turbulence  data  clearly  indicate  that  the  probe  volume  was 
traversed  from  a  point  upstream  of  the  shock  to  a  point 
downstream  of  the  shock.  The  measured  velocity  decreased  from 
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Figure  27 .  Recalculated  Data  for  Boundary  Layer  Survey 
031693a,  Wall  to  Free-Stream 
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Figure  28.  Recalculated 
Wall  to  Free-Stream 
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Figure  32.  Velocity  and  Turbulence  Intensity  Plots,  Shock 
Survey  031093b 
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a  known  free-stream  supersonic  value  to  a  known  free-stream 
subsonic  value  as  the  probe  volume  was  traversed  across  the 
normal  shock  position.  Also,  turbulence  intensity  increased 
from  free-stream  conditions  as  the  probe  volume  crossed  the 
shock  then  decreased  to  free-stream  values  downstream  of  the 
shock.  The  dramatic  increase  in  turbulence  intensity  (above 
10%)  indicated  a  region  of  flow  unsteadiness.  Between  0.8  and 
1.2  inches  on  survey  030393a  and  between  1.0  and  1.4  inches  on 
survey  031093b  there  was  a  region  of  decreased  mean  particle 
velocity  and  double  peaked  histograms. 

This  region  was  interpreted  as  the  region  of  unsteadiness 
of  the  shock.  The  schlieren  visualization  for  both  runs,  as 
shown  in  Figures  33  and  34,  supported  this  conclusion.  The 
random  photographs  taken  immediately  prior  to  each  LDV  survey 
indicated  the  shock  to  be  fluctuating  between  0.75  and  1.25 
inches  during  030393a  and  between  1.0  and  1.4  inches  during 
031093b.  This  corresponded  very  closely  to  the  information 
taken  from  the  plot  of  the  LDV  data. 

In  the  region  where  the  probe  volume  was  traversed  through 
the  shock,  the  LDV  data  became  difficult  to  interpret  due  to 
the  presence  of  double  peaked  histograms.  The  points  in  this 
region  on  the  plots  in  Figures  31  and  32  represent  the  mea*' 
velocity  taken  from  the  double  histogram.  The  appearance  of 
the  double  histogram  in  this  area  was  concluded  to  be  a  result 
of  the  shock  movement  or  unsteadiness.  As  the  normal  shock 
moved  back  and  forth  across  the  stationary  probe  volume,  the 
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LDV  system  would  alternately  measure  supersonic  and  subsonic 
particle  velocities.  The  relative  strength  of  each  peak  of  the 
histogram  would  logically  be  related  to  the  the  length  of  time 
the  probe  volume  spent  in  either  the  subsonic  or  supersonic 
region  of  the  flow.  This  conclusion  is  supported  by 
experimental  data  obtained  by  Strazisar  [Ref. 22]  and  Chriss 
et.al.  [Ref  23]  in  which  double  histograms  were  observed 
during  LDV  surveys  through  a  normal  shock. 

Figure  35  shows  an  additional  plot  of  the  LDV  data  taken 
in  031093b  in  which  an  attempt  was  made  to  further 
characterize  the  velocity  distribution  across  the  region  of 
unsteadiness.  In  each  survey,  the  points  where  double 
histograms  appeared  were  separately  analyzed  by  splitting  the 
double  histograms  into  a  low  and  a  high  peak.  By  editing  each 
double  histogram,  a  mean  value  for  each  peak  was  determined. 
The  amount  of  variance  of  each  high  and  low  mean  was 
calculated,  and  this  was  an  indication  of  the  turbulence 
intensity.  The  separation  of  the  histograms  was  a  measure  of 
the  flow  unsteadiness.  Figure  35  shows  the  high  and  low  means 
and  the  amount  of  variance  associated  with  each  point  inside 
the  region  of  unsteadiness.  Ideally,  a  more  rational 
mathematical  method  of  splitting  the  histograms  by  using  the 
statistical  information  is  required  in  order  to  obtain  a  full 
understanding  of  the  velocity  information  obtained  in  the 
double  peaked  histograms,  as  apposed  to  the  manual  editing 
technique  used  here. 
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One  very  repeatable  phenomenon  associated  with  the  normal 
shock,  surveys  was  the  region  just  downstream  of  the  shock 
movement  where  the  velocity  bottomed  out  before  rising  to  a 
free- stream  value.  This  can  be  seen  in  Figures  31  and  32.  This 
behavoir  agrees  with  That  reported  by  Chriss  et.al.  [Ref.  23] 
who  hypothesized  that  it  was  due  to  the  acceleration  of  the 
free  stream  flow  as  a  result  of  the  significant  growth  of  the 
boundary  layers  on  all  four  walls. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  analysis  of  the  data  in  this  study  leads  to  a  number 
of  conclusions  about  LDV  measurements  in  a  transonic  flow  and 
about  particle  dynamics  across  a  normal  shock: 


*  LDV  measurements  made  upstream  and  downstream 
of  the  shock  were  very  successful  and  provided 
accurate  and  repeatable  velocity  magnitudes  in 
backscatter  mode  in  a  transonic  flow. 

*  Seeding  at  the  contraction  with  1  fj.m  sized  PSL 
particles  proved  adequate.  Data  rates  were  more 
than  sufficient  in  all  phases  of  the  experiment. 

*  The  design  of  the  supersonic  wind  tunnel  nozzle  blocks 
was  flawed  in  that  insufficient  expansion  was  provided 
to  compensate  for  boundary  layer  effects.  The  tunnel 
must  presently  be  run  -with  back  pressure  relief  in 
order  to  allow  movement  of  the  starting  normal  shock 
downstream  of  the  test  section. 

*  The  design  of  the  plenum,  specifically  the  absence  of 
any  screens  or  flow  straightening  devices,  caused 
turbulence  which  might  have  contributed  to  shock 
unsteadiness . 

*  Shock  position  and  movement  could  be  accurately 
determined  with  schlieren  visualization  and  LDV 
measurements.  Spark  source  schlieren  photography 
proved  to  be  successful  as  a  means  of  recording  the 
extent  of  shock  movement.  Data  obtained  from  LDV 
measurements  also  proved  to  be  accurate  in  the 
determination  of  shock  location. 
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*  There  was  an  apparent  region  at  the  downstream  edge 
of  the  normal  shock  movement  where  particle  velocity 
reached  a  minimum  before  accelerating  to  free  stream 
velocity.  This  phenomenon  was  highly  repeatable  and 
may  be  a  function  of  the  downstream  boundary  layer 
growth. 

*  Instability  of  the  normal  shock  caused  an  area  of 
double  peaked  histograms  which  can  possibly  be 
statistically  analyzed  to  obtain  particle  dynamic 
information  through  the  normal  shock. 


B .  RECOMMENDATIONS 

LDV  is  a  powerful  tool  with  the  ability  to  provide 
accurate  velocity  measurements  under  specific  conditions.  The 
accuracy  of  the  data  does,  however,  depend  a  great  deal  on  the 
ability  to  fully  understand  and  correctly  model  particle 
dynamics.  More  experimentation  and  improvements  in  the  current 
test  facilities  and  apparatus  are  required  to  continue  LDV 
work  in  the  supersonic  tunnel.  The  following  recommendations 
are  made : 


*  The  supply  air  pressure  must  be  increased  at  the 
facility  to  provide  longer  run  times.  Certifying 
additionaj.  storage  tank  would  also  increase  run 
time  by  increasing  the  storage  volume.  Longer  run 
times  are  the  key  to  obtaining  more  detailed  surveys. 

*  The  control  and  monitoring  of  the  supersonic  tunnel 
(supply  valve,  plenum  pressure  gauge)  should  be 
incorporated  into  a  master  control  panel  which  would 
allow  better  management  of  the  tunnel  operation. 

*  Reduction  of  turbulence  can  be  accomplished  by  placing 
screens  or  other  flow  straightening  devices  in  the 
plenum. 
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*  The  HP  data  acquisition  system  can  be  utilized  further 
by  incorporating  thermal  compensation  cards  and  the 
ability  to  measure  real-time  values  of  plenum  and  test 
section  temperature. 


*  The  current  traverse  table  for  the  LDV  optics  is 
unable  to  fully  traverse  the  test  section  and  must 
be  operated  manually.  An  automatic  traverse  system 
with  position  indication  would  allow  more  accurate 
surveys . 
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APPENDIX  A.  WIND  TUNNEL  DATA  ACQUISITION  SYSTEM 


A.  DESCRIPTION  OP  SYSTEM 

The  data  acquisition  system  for  the  supersonic  wind  tunnel 
consists  of  the  following  items:  2  pressure  transducers,  HP 
9000  series  300  computer  data  acquisition/reduction  system,  HP 
3455A  digital  voltmeter,  HP  3497A  data  acquisition/control 
unit,  Lockheed  MT-1758  test  panel,  and  "SPEED_5"  software.  The 
purpose  of  the  data  acquisition  is  to  measure  total  plenum 
pressure  (P2)  and  static  test  section  pressure  (Px) .  This 
information  is  then  used  to  calculate  test  section  Mach  number 
upstream  of  the  normal  shock  based  on  the  following  isentropic 
flow  equation: 

El  =  [i  +  y  -  1m2]^~t 

Pl  2 

Test  section  stagnation  temperature  (  measured  in  the  plenum 
and  input  by  the  user)  and  Mach  number  are  then  used  to 
calculate  velocity  upstream  of  the  shock.  Finally,  by  dividing 
velocity  by  the  known  fringe  spacing,  a  Doppler  frequency  is 
calculated. 

Once  the  system  is  calibrated  and  started,  all 
measurements  and  calculations  are  automatically  and 
continuously  performed  throughout  the  entire  wind  tunnel  run. 
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A  continuous  display  of  P2,  Px,  P2/Plf  Mach  number,  and 
frequency  are  presented  on  the  computer  CRT  and  sent  to  the 
printer.  The  presented  data  can  be  utilized  during  the  run  to 
monitor  tunnel  conditions  and  properly  set  frequency  filters 
on  the  processor.  The  data  is  also  useful  as  a  means  of 
recording  a  time  history  of  tunnel  conditions  for  each  run. 

The  "SPEED_5"  program  was  developed  from  the  "READ_ZOC" 
utility  program  given  in  Reference  24.  A  small  portion  of  the 
original  program  was  adapted  for  use  with  the  supersonic  wind 
tunnel  acquisition  system. 

B.  OPERATING  THE  DATA  ACQUISITION  SYSTEM 

1.  Start-up 

*  Turn  on  the  HP9000,  HP3455A,  and  HP3497A. 

*  Verify  supply  air  is  connected  to  MT-1758  and 
calibration  pressure  is  set,  using  the  high 
pressure  regulator  knob,  to  50.9  inches  Hg 
(25  psi) . 

*  Type  SHIFT-RESET  to  enter  HP  BASIC. 

*  LOAD  "SPEED_5 " 

*  RUN  "SPEED_5 " 

2.  Calibration 

*  Type  0  to  continue  program  after  the  introduction. 

*  ENTER  the  atmospheric  pressure  in  psi. 

*  ENTER  the  test  section  static  temperature 
(calculated  from  previous  runs  where  plenum 
stagnation  temperature  was  directly  measured) . 

*  ENTER  0  to  calibrate  P^  Open  Px  to  atmosphere  and 
zero  reading  on  HP3455A  then  apply  calibration 


69 


pressure  and  adjust  to  12.5  on  HP3455A. 

*  ENTER  4  and  repeat  above  process  to  calibrate  P2. 

*  Type  an  out  of  range  value  (any  number  other  than 
0  through  9)  to  run  the  program. 

*  The  program  will  run  continuously  until  SHIFT-RESET 

is  typed. 


A  listing  of  the  "SPEED_5"  program  can  be  seen  in  Figure 

A1 . 
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APPENDIX  B.  TSI  "FIND"  SOFTWARE 


TSI  Flow  Information  Display  (FIND)  software  was  used  as 
the  data  acquisition  interface  for  the  LDV  system.  The 
software  is  designed  to  acquire,  analyze,  and  reduce  LDV  raw 
data.  FIND  also  allows  the  user  to  establish  and  set  up 
hardware  parameters.  Three  subprograms  (  data  acquisition, 
statistical  data  analysis,  and  flow  field  plot)  were 
extensively  used  during  the  course  of  this  study.  Figures  B1 
through  B8  show  a  typical  sequence  of  displays  within  the  FIND 
software  used  to  set  up  and  run  an  LDV  survey. 

The  main  menu  is  shown  in  Figure  B1 .  From  the  main  menu, 
"A"  is  selected  to  access  the  data  acquisition  program  which 
is  shown  in  Figure  B2 .  From  Figure  B2 ,  "I",  "P",  and  "0"  are 
selected  to  set  the  I/O  port  and  processor  selections, 
processor  set  up,  and  optics  configurations.  Typical 
selections  and  settings  are  shown  in  Figures  B3,  B4,  and  B5 . 
To  complete  the  set  up,  "F"  is  selected  from  the  main  menu  to 
enter  the  data  files  management  screen  as  shown  in  Figure  B6 . 
After  acquiring  the  raw  data  from  an  LDV  survey,  the 
Statistical  Analysis  and  Flow  Field  Plot  programs  are  selected 
as  shown  in  Figures  B7  and  B8  to  analyze  and  reduce  the  data. 
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Figure  B4.  "FIND"  Processor  Set  Up  Screen 
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Figure  B5 .  "FIND"  Optics  Configuration  Screen 


Figure  B6 .  "FIND"  Data  Files  Management  Screen 
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Figure  B7 .  "FIND"  Statistical  Analysis  Program  Menu 


Figure  B8.  "FIND"  Flow  Field  Plot  Program  Menu 


77 


APPENDIX  C.  SUPERSONIC  FREE- JET 


A  supersonic  free  jet  was  designed  and  constructed  during 
the  course  of  this  study.  The  intent  was  to  have  available  a 
means  of  producing  supersonic  flow  in  the  event  that  the 
supersonic  blow  down  tunnel  proved  inoperable  due  to  design 
errors  or  unusable  due  to  limited  run  times.  The  free  jet,  as 
shown  in  Figure  Cl,  is  capable  of  essentially  continuous 
operation  with  supply  pressure  at  150  psi  and  has  the 
advantage  of  direct  optical  adcess. 

Supply  air  is  routed  directly  to  a  regulator  which  serves 
as  the  control  valve  for  the  free  jet.  The  regulator  connects 
directly  with  an  4  inch  diameter  straight  pipe  that  serves  as 
the  plenum.  The  regulator  is  also  supplied  with  back  pressure 
from  the  plenum  which  it  uses  to  maintain  the  pressure  within 
the  plenum.  A  solid  aluminum  nozzle  is  attached  to  the  end  of 
the  plenum.  Air  exits  the  nozzle  and  is  captured  a  few  feet 
away  by  conical  ducting  which  routes  the  discharge  to  the 
exterior  of  the  building.  A  pressure  gauge  is  attached  to  the 
plenum  to  monitor  operating  conditions  and  additional  ports 
are  available  to  measure  temperature  and  allow  seeding. 

The  aluminum  nozzle  is  an  ASME  low  0  series,  long- radius 
flow  nozzle  [Ref  25] .  The  coefficient  0  represents  the  ratio 
of  nozzle  exit  diameter  to  inlet  diameter  where  low  0  indicates 
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Figure  Cl.  Supersonic  Free 


a  ratio  below  0.5.  The  free- jet  nozzle  exit  was  chosen  as 
1.25"  giving  it  a  diameter  ratio  of  about  0.31.  The  remaining 
dimensions  of  the  nozzle,  including  the  elliptical  approach  to 
the  throat,  are  determined  by  the  choice  of  exit  diameter. 
Throat  length  is  defined  as  0.6  times  the  exit  diameter.  The 
minor  axis  of  the  elliptical  approach  is  defined  as  0.66  times 
the  exit  diameter  and  the  major  axis  is  defined  as  equal  to 
the  exit  diameter.  The  design  of  the  nozzle  is  further 
illustrated  in  Figure  C2 .  Although  the  free  jet  was  not 
utilized  during  this  study,  it  is  a  useful  addition  to  the 
laboratory  in  that  it  will  allow  further  LDV  studies  to  be 
performed  under  continuous,  predictable,  and  highly 
controllable  conditions. 
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Figure  C2 .  Supersonic  Free- Jet  Nozzle  Design,  From 
Ref.  24 
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APPENDIX  D.  TABULATED  RESULTS 


The  following  is  a  compilation  of  the  results  from  the 
present  study.  The  data  are  presented  in  the  following  format: 


*  BOUNDARY  LAYER  SURVEY  RAW  DATA 

*  BOUNDARY  LAYER  SURVEY  PLOTS  OF  MEAN  VELOCITY  AND 
TURBULENCE  VERSUS  POSITION 

*  SHOCK  SURVEY  RAW  DATA 

*  SHOCK  SURVEY  PLOTS  OF  MEAN  VELOCITY  AND  TURBULENCE 
INTENSITY  VERSUS  POSITION 
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